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We discuss the dynamics of a Bose-Einstein condensate in a double- well trap subject to phase noise 
and particle loss. The phase coherence of a weakly-interacting condensate, experimentally measured 
via the contrast in an interference experiment, as well as the response to an external driving become 
maximal for a finite value of the dissipation rate matching the intrinsic time scales of the system. 
This can be understood as a stochastic resonance of the many-particle system. Even stronger effects 
are observed when dissipation acts in concurrence with strong inter-particle interactions, restoring 
the purity of the condensate almost completely and increasing the phase coherence significantly. Our 
theoretical results are backed by Monte Carlo simulations, which show a good qualitative agreement 
and provide a microscopic explanation for the observed stochastic resonance effect. 

PACS numbers: 03.75.Lm, 03.75.Gg, 03.65.Yz 



I. INTRODUCTION 

Stochastic resonance (SR) is a strongly surprising, 
yet very general effect in nonlinear dynamical systems. 
Against our naive understanding, the response of a sys- 
tem to an external driving can be facilitated if an appro- 
priate amount of noise is added. In fact, the maximum 
of the response - the stochastic resonance - is found if 
the timescale of the noise matches an intrinsic time scale 
of the system. The effect was first described for strongly 
damped classical systems such as the overdamped parti- 
cle in a driven double well trap. In this case the noise 
is strong enough to induce the transition between the 
wells, whereas it is still weak enough not to randomize 
the dynamics completely. The particle will then hop to 
and fro almost deterministically if the average transition 
time between the wells due to the noise equals half of 
the driving period By now, a stochastic resonance 
has been shown in a variety of systems, an overview can 
be found in the review articles [1, H, III HI • In addition 
to numerous examples in classical dynamics, stochastic 
resonance has also been found in a variety of quantum 
systems (see, e.g. , d S, 0, i i, M ED] ) ■ 

Recently, there has been an increased interest in the ef- 
fects of dissipation and the possibilities to control these in 
interacting many-body quantum systems. For instance, 
the entanglement in a spin chain assumes an SR-like max- 
imum for a finite amount of thermal noise fl2| . Methods 
to attenuate phase noise for an open two-mode BEC were 
discussed in and the effects of particle loss on the 
spin squeezing of such a system were analyzed in [l4| . 
Furthermore, it has been shown that dissipative processes 
can be tailored to prepare arbitrary pure states for quan- 
tum computation and strongly correlated states of ultra- 
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cold atoms [la, Il6l or to implement a universal set of 
quantum gates ,17:]. Actually, a recent experiment has 
even proven that strong inelastic collisions may inhibit 
particle losses and induce strong correlations in a quasi 
one-dimensional gas of ultracold atoms [H, [l|| . 

In the present paper, we investigate the effects of noise 
and dissipation for a Bose-Einstein condensate (BEC) in 
a double-well trap. The essential idea has been intro- 
duced in a recent letter [2(|, and here we extend the 
discussion to a detailed analysis of the predicted SR- 
phenomenon. The setup under consideration has been 
experimentall y re alized by several groups only in the last 
few years [H [H [H, H, Hi, HI]. Ultracold atoms in op- 
tical potentials have the enormous advantage that they 
allow to observe the quantum dynamics of an interacting 
many-particle system in situ. Thus they serve as excel- 
lent model systems, bringing together aspects of nonlin- 
ear dynamics, solid-state physics and the theory of open 
quantum systems. Here we show that the coherence of 
the two condensate modes assumes a maximum in the 
fashion of the stochastic resonance effect for a finite dissi- 
pation rate, which matches the time scales of the intrinsic 
dynamics. In this case the particle loss is strong enough 
to significantly increase the condensate purity, whereas it 
is still weak enough not to dominate the dynamics com- 
pletely. Similarly the response to an external driving is 
increased if a proper amount of dissipation is present. 
Even more remarkable results are found when dissipa- 
tion acts in concurrence with strong inter-particle inter- 
actions, leading to an almost complete revival of the pu- 
rity of the BEC. These effects are of considerable strength 
for realistic parameters and thus should be readily ob- 
servable in ongoing experiments. 

This paper is organized as follows: First, we introduce 
the theoretical description of the open two-mode Bose- 
Hubbard system. We discuss the main sources of noise 
and dissipation and derive the corresponding mean-field 
approximation of the many-particle system. The result- 
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ing dynamics for weak inter-particle interactions is an- 
alyzed in Sec. IIIII It is shown that the phase contrast 
between the two modes assumes an SR-like maximum if 
the time scales of tunneling and dissipation are matched. 
This result is explained within the mean-field approxima- 
tion as well as for the underlying many-particle quantum 
dynamics with Monte Carlo simulations backing up the 
approximative results. The response of the open system 
to an external driving is discussed in Sec. |W] The am- 
plitude of the forced oscillation also shows a pronounced 
stochastic resonance effect. Sec. IVlthen investigates the 
case of a strongly interacting BEC, which is a problem of 
both, fundamental theoretical interest as well as high ex- 
perimental relevance. The interplay between interactions 
and dissipation can restore the purity of the condensate 
almost completely and significantly increase the phase co- 
herence in comparison with situations where one of the 
two is weak or missing. This counter-intuitive effect is 
robust and can be explained by the appearance of novel 
nonlinear eigenstates. 

II. NOISE AND DISSIPATION IN A TRAPPED 

BEC 

The basic setup under consideration is depicted in 
Fig. [T] Ultracold atoms are confined in a double-well 
trap that can be realized e.g. by superimposing an op- 
tical lattice with an optical dipole trap [2l|, |22|, [23| , in a 
bichromatic optical lattice [24], (2|| , or on an atom chip 
[26| . We consider the case where only one mode in each 
trap is significantly populated, whereas all higher modes 
contribute to the heat bath (see below). The unitary dy- 
namics of the atoms is then described by the two-mode 
Bose-Hubbard Hamiltonian [27|, H, [H, [3(J 

H = -J (a\a 2 + a\a^j + e 2 h 2 + ei<h\ 

+-(ni(ni-l) + n 2 (n 2 -l)) ) (1) 

where a 7 - and a\ are the bosonic annihilation and creation 

operators in mode j and fij = a^cij is the corresponding 
number operator. Furthermore, J denotes the tunnel- 
ing matrix element between the wells, U the interaction 
strength and Cj the on-site energy of the jth well. We set 
H = 1, thus measuring all energies in frequency units. 

In order to clarify the algebraic structure of the model 
and to analyze the dynamics in the Bloch representation 
we introduce the collective operators 

L x = - (a\a 2 + a\a^j 

L y = \ (a!«2 - ( 2 ) 

L z = i [a\a 2 - a\a^ , 

which form an angular momentum algebra su(2) with 
quantum number i = N/2 [13, [H, H§, H HI, HI , where 




FIG. 1: (Color online) The open double-well trap considered 
in the present paper. 

N is the actual particle number. The Hamiltonian ((T|) 
then can be rewritten as 

H = -2JL X + 2eL z + UL 2 Z (3) 

up to terms only depending on the total number of atoms. 
Here, e = e 2 — t\ denotes the difference of the on-site 
energies of the two wells. 

A model for noise and dissipation in a deep trapping 
potential has been derived by An glin [33| and later ex- 
tended by Ruostekoski and Walls [34[ to the case of two 
weakly coupled modes. The dynamics is then given by 
the master equation 

'p = -i[H, p}-yY, (tfp + P&j ~ 2f w) ( 4 ) 

.7=1,2 

-y E (c} ± c 3±P + pc} ± c 3± -2c ]±P c} ± ) 

with the Lindblad operators 

Cj+ = and 

C s - = e^^-^+^aj, (5) 

describing growth and depletion of the condensate. 

Let us briefly discuss the effects of the noise and dis- 
sipation terms. The second term ~ j p in Eqn. (|4|) 
describes phase noise due to elastic collisions with the 
background gas atoms. It is usually the dominating 
contribution, effectively heating the system, but leav- 
ing the total particle number invariant. If only phase 
noise is present, the system relaxes to an equilibrium 
state where all coherences are lost and all Dicke states 
I rii, N — ill) ~ a\ ni al N ~ ni |0, 0) are equally populated 

(ni,N- mlpln'^N- n[) = N 5 Wl , n /, (6) 

as long as J ^ [H, HH . This corresponds to a thermal 
state of infinite temperature with (L) = 0. The remain- 
ing terms ~ "f a in the master equation Q describe am- 
plitude noise, i.e. the growth and depletion of the con- 
densate due to inelastic collisions with the background 
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gas. In contrast to phase noise, amplitude noise heats 
and cools the system. If both amplitude and phase noise 
are present, the system will relax to the proper thermal 
state with a density operator p oc exp(—f3(H — ph)) [33| . 

In current experiments amplitude noise and dissipation 
is usually extremely weak in comparison to phase noise 
[34| . if it is not introduced artificially as for example by 
forced evaporative cooling during the preparation of the 
BEC. For example, phase noise damps Josephson oscilla- 
tions on a timescale of a few hundred milliseconds in the 
experiments, while less than 10% of the atoms are lost 
during a 30 s experiment [U [22|, [23| . This is much too 
weak to produce the effects discussed in the present pa- 
per. In contrast, a strong and tunable source of dissipa- 
tion can be implemented artificially by shining a resonant 
laser beam onto the trap, that removes atoms with the 
site-dependent rates 7<y from the two wells j = 1,2. In 
magnetic trapping potentials, this can also be achieved 
by a forced rf-transition to an untrapped magnetic sub- 
state 37] . Non-trivial effects of dissipation such as the 



stochastic resonance discussed below require strongly bi- 
ased loss rates, i.e. j a i ^ 7a2- For a laser beam focused 
on one of the wells an asymmetry of f a = (702 - 7ai)/(7o2 
+ 7a 1) = 0.5 should be feasible. The master equation de- 
scription of both, phase noise and particle loss, is well es- 
tablished (3II ] and routinely used in the context of photon 
fields. In the following we will thus consider the dynamics 
generated by the master equation 

p = -i{H 7 p}~^ (n 2 jP + Pn? -2hjphj) 
3=1,2 



9 la i (a]ajp + pa]aj - 2a JP af\ . (7) 



J=l,2 

The macroscopic dynamics of the atomic cloud is to 
a very good approximation [2^, |3(| IH[ described by a 
mean-field approximation, considering only the expec- 
tation values Sj(t) = 2ti(Ljp(t)) of the angular mo- 
mentum operators ([2]) and the particle number n(t) = 
tr((rii -I- fi2)p(t)). The evolution equations for the Bloch 
vector s = (s x ,s y ,s z ) are then calculated starting from 
the Master equation via Sj = tr(Ljp) with the exact re- 
sult (cf. H) 



2Js z + 2es x + U(s x s z + A 



= -2Js y - T t 1 s z - T t 1 f a n 1 



-T, 



v 

-1„ rr-l 



Ti fa,S z , 



(8) 



where we have defined the transversal and longitudinal 
damping times by 



Tf 1 = ( 7 ai +7a 2 )/2 and T 2 ~ l = lp + T{ 



(9) 



These equations of motion resemble the celebrated Bloch 
equations in nuclear magnetic resonance [39T |40T |. with 
some subtle but nevertheless important differences. The 



longitudinal relaxation is now associated with particle 
loss and, more important, the dynamics is substant ially 
altered by the U -dependent interaction term [2l|, [53, l28f . 

The exact equations of motion © still include the co- 
variances 



Aj/c — (LjLk + Lf~Lj) — 2(Lj)(Lk}- 



(10) 



The celebrated mean-field description is now obtained by 
approximating the second order moment s by products of 
expectation values such that A jfc w [U [H, [29|, |30T |. 
This truncation is valid in the macroscopic limit of large 
particle numbers, since the covariances vanish as 1/n if 
the many particle quantum state is close to a pure BEC. 

In the following, we will show that a finite amount of 
dissipation induces a maximum of the coherence which 
can be understood as an stochastic resonance effect. In 
this discussion we have to distinguish between two dif- 
ferent kinds of coherence, which will both be considered 
in the following. First of all we consider the phase co- 
herence between the two wells, which is measured by the 
av erag e contrast in interference experiments as described 
in [11 [Hil] and given by 



«(*) 



2|(&la 2 )| _ y/ Sx (t) 2 + s y {ty 

(fix + n 2 ) 



n(t) 



(11) 



Secondly, we will analyze how close the many-particle 
quantum state is to a pure Bose-Einstein condensate. 
This property is quantified by the purity 



p = 2tr(p2 cd )-l 



(12) 



of the reduced single-particle density matrix 29, 30, 32, 

El 



Prod 



1 

N 



(a[ai) (a|a 2 ) 
(ajai) (4a 2 ) 



(13) 



One can easily show that the purity is related to the Bloch 
vector sbyp=|s| 2 /n 2 . A pure BEC, corresponding to 
a product state, is, of course, characterized by p = 1. 



III. DISSIPATION INDUCED COHERENCE IN 
A WEAKLY-INTERACTING BEC 

In this section, we show that a proper amount of dissi- 
pation can indeed increase the phase coherence (fTTj) of a 
two- mode BEC similar to the stochastic resonance effect. 
For simplicity, we start with the linear case U = 0, where 
the mean-field equations of motion for the expectation 
values (|5J) are exact. The linear equations resemble the 
Bloch equations for driven nuclear spins in the rotating 
wave approximation [40| , which are known to show a pro- 
nounced stochastic resonance effect [39| : The amplitude 
of forced oscillations of the spins given by s y assumes a 
maximum for a finite value of the relaxation rates T^ 1 
and T 2 _1 , provided these are coupled. For the two-mode 
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BEC considered here this is automatically the case as 
given by Eqn. @. Thus we also expect a maximum of 
the steady state value of the phase coherence (fTT|) for a 
finite value of T-f" . 

Let us now determine the steady state value of the 
contrast (fTTj) which quantifies the phase coherence of the 
two wells, as a function of the system parameters and 
the relaxation rates. Obviously, the only steady state in 
the strict sense is given by s = and n = 0, correspond- 
ing to a completely empty trap. However, the system 
rapidly relaxes to a quasi-steady state where the internal 
dynamics is completely frozen out and all components of 
the Bloch vector and the particle number decay at the 
same rate: 



(a) 



s(t) = s e~ 



n(t) = noe 



(14) 



Substituting this ansatz into the equations of motion ([8]) , 
the quasi-steady state is determined by the eigenvalue 
equation 



M 



with the matrix 



M 



(T 2 
-2e 


V o 



/s x o\ 




/Sx0\ 


SyQ 


- 


SyO 


SzO , 




1 Sz0 J 


\n J 




\n J 


1 2e 








(15) 



T, 



i 
2,1 



-21 
TV 1 






f a Tr 



(16) 



faT,- 1 Tf 1 J 



which is readily solved numerically. 

Fig. [2] shows the resulting values of the contrast a as 
a function of the dissipation rate T^ 1 and the tunneling 
rate J for U = e = and 7 p = 5s _1 . For a fixed value 
of one of the parameters, say J, one observes a typical 
SR-like maximum of the contrast for a finite value of the 
dissipation rate X/T\ as shown in part (b) of the figure. 
In particular, the contrast is maximal if the time scales of 
the tunneling and the dissipation are matched according 
to 



(17) 



Furthermore, the contrast a(J) shows a similar maxi- 
mum for a finite value of the tunneling rate J when the 
dissipation rate is fixed as shown in Fig. [5] (c). Con- 
trary to our intuition this shows that an increase of the 
coupling of two modes can indeed reduce their phase co- 
herence. 

In the special case e = 0, illustrated in Fig. [21 one 
can solve the eigenvalue equation (fTS"]) exactly. In this 
case one has s x = and the contrast a is related to the 
eigenvalue k by 



2J(Tf 1 - k) 
faT^iT^-n) 



(18) 




1^ [S-] 



FIG. 2: (Color online) Contrast a of the quasi-steady state 
(|14p as a function of the tunneling rate J and the dissipation 
rate 1/Ti (a) for 7 P = 5 s _1 and U = e = and (b) for a fixed 
value of the tunneling rate J = 2 s" 1 and (c) a fixed value of 
the dissipation rate 1/Ti = 2s _1 . The dash-dotted red lines 
represent the approximations (|19[) for small and large values 
of J, respectively. 



Evaluating the roots of the characteristic polynomial to 
determine K leads to an algebraic equation of third order 
which can be solved analytically. The resulting expres- 
sions are quite lengthy, but the limits for small and large 
values of the tunneling rate are readily obtained as 



2 J 



Tf L -(l-f a )T~ L 

faT' 1 

rv sa — — 1 — 

" 2,7 



for J <C T-f 
for J >Tf 



(19) 



These approximations are plotted as dashed red lines in 
Fig. [5] (c). Their intersection given by (fTT|) gives a very 
good approximation for the position of the SR-like max- 
imum of the contrast a{J). 

An important experimental issue is the question 
whether the quasi-steady state is reached fast enough, 
such that the typical SR-like curve of the contrast as 
shown in Fig. [21 can be observed while still enough atoms 
are left in the trap. To answer this question, we integrate 
the equations of motion (jSJ) starting from a pure BEC 
with s(0)/n(0) = (V3/2, 0,1/2) and n(0) = 100 parti- 
cles. Fig. [3] (a) shows the relaxation of the contrast for 
J = 4s _1 and T\ = Is. The steady state value is nearly 
reached after t — 1 s, when still 40% of the atoms are 
left in the trap. Fig. [3] (b) shows the development of the 
contrast a(J) in time. It is observed that the character- 
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(a) 



FIG. 3: (Color online) Relaxation to the quasi-steady state 
for 7 P = 5s" 1 , Tf 1 = Is" 1 , e = 10s" 1 and U = 0. (a) 
Relaxation of the contrast a(t) for J = 4 s" 1 . (b) Decay of 
the particle number n(t) for J = 4 s" 1 . (c) Development of 
the SR-maximum of the contrast a(J). 



istic SR-like maximum is already well developed after 1 
second, where roughly half of the atoms are lost. Thus 
we conclude that the SR-like maximum of the contrast 
should be observable in ongoing experiments. 

The stochastic resonance effect introduced above is ro- 
bust and generally not altered by changes of the system 
parameters or in the presence of weak inter-particle inter- 
actions. For instance, a change of the bias e of the on-site 
energies of the two wells preserves the general shape of 
a(l/Ti, J) shown in Fig. [21 and especially the existence 
of a pronounced SR-like maximum. At most, the func- 
tion cv(l/Ti, J) is stretched, shifting the position of the 
SR-like maximum. This shift is illustrated in Fig. 2] (a) 
where we have plotted the contrast as a function of J 
for the dissipation rate = 2 s" 1 and different values 
of e. Thus, this effect provides a useful tool to shift the 
maximum to values of J, which are easier accessible in 
ongoing experiments. 

Similarly, the position of the maximum of the coher- 
ence a (J) is shifted in the presence of weak inter-particle 
interactions. An interacting BEC will usually not show 
a simple exponential decay of the form ([14]) because the 
instantaneous decay rate depends on the effective inter- 
action strength Un(t), which also decreases [42l. l43l |44|. 
However, the discussion of quasi-steady states and in- 
stantaneous decay rates is still useful if the decay is weak. 
In this case the system can follow the quasi-steady states 
adiabatically and the decay of the population is given by 



dn{t) 

dt 
ds(t) 

dt 



= -K(n(t))n(t) 
= -/e(n(t))s(i) 



and 



(20) 



in good approximation. Substituting this ansatz into the 
equations of motion {5| yields four coupled nonlinear al- 
gebraic equations, which can be disentangled with a little 




J Is" 1 ] 



FIG. 4: Steady state values of the contrast a as a function 
of the tunneling rate, for (7 = and different values of the 
energy bias e (a), and as a function of the effective interaction 
strength g — Un for e = (b). The remaining parameters are 
7 P = 5s" 1 and Tf 1 = 2s" 1 . 



algebra. For a given number of particles n, the instan- 
taneous decay rate k is obtained by solving the fourth 
order algebraic equation 



U z ftT^{n-T~'){n-T 2 



(21) 



The Bloch vector for the corresponding quasi-steady 
state is then given by 



SyQ 



SzO 



Tr 1 (k - T{ 



1\2 



-p a Tr 



K — T 



rr 1 ) 2 



2J/ Q 2 Tf 2 



■Un 2 



k - Tf 1 



f a Tr 



(22) 



The fourth order equation (|2ip yields four solutions for 
the decay rate k. Discarding unphysical values, one finds 
either one or three quasi-steady states. This appearance 
of novel nonlinear stationary states has been discussed in 
detail in the context of nonlinear Landau-Zener tunneling 
[H, HI, S3, Ell and nonlinear transport [49l [50l] . 

The resulting contrast a (J) in a quasi-steady state is 
shown in Fig. 0] (b) for different values of the effective 
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FIG. 5: (Color online) (a) Histogram of the probabilities to 
measure the relative phase <j> and the relative population im- 
balance s z in a single experimental run after t = 1.5 s obtained 
from a MCWF simulation of the many-body dynamics. The 
initial state was chosen as a pure BEC (i.e. a product state) 
with s z = n/2 and n(0) = 100 particles and the remaining pa- 
rameters are 7 P = 5s , T\ = 0.5s, e = 10s _1 , U — 0.1s -1 . 
(b) Average contrast a = ^/s| + s|/n (solid black line) af- 
ter t = 1.5 s compared to \/s| + sf/|s| and \s\/n (dashed red 
lines). 



thereby given by 

P(s z ) = tr(|s z ) (s z \p) and 
P(4>) = tr(|0> <0| p), 

where the L z eigenstates 

\s z ) = | n/2 - s z ,n/2 + s z ) with 
s 2 = -n/2,-7i/2+l,...,n/2 

and the phase eigenstates 

+n/2 

10):=^== J] e*"| s ,) with 

1 n 

4> = 0, 2tt , . . . , 2tt 

K ' n+1' n+ 1 



(23) 



(24) 



(25) 



each form a complete basis. 

Part (a) of Fig.[5]shows a histogram of the probabilities 
to observe the relative population imbalance s z /n and the 
relative phase <p in a single experimental run for three 
different values of the tunneling rate J after the system 
has relaxed to the quasi-steady state. With increasing J, 
the atoms are distributed more equally between the two 
wells so that the single shot contrast increases. Within 
the mean-field description this is reflected by an increase 

of * s% + s%/\s\ at the expense of |s z |/|s| (cf. part (b) of 



the figure). However, this effect also makes the system 
more vulnerable to phase noise so that the relative phase 
of the two modes becomes more and more random and 
|s|/n decreases. The average contrast (jTTj) then assumes 
a maximum for intermediate values of J as shown in part 
(b) of the figure. 



IV. STOCHASTIC RESONANCE OF A DRIVEN 

BEC 



interaction constant g = Un. One observes that the po- 
sition of the SR-like maximum of the contrast is shifted 
to larger values of the tunneling rate, while the height re- 
mains unchanged. Furthermore the shape of the stochas- 
tic resonance curve a (J) is altered, becoming flatter for 
J < >/max and steeper for J > J max . For even larger 
values of the interaction constant Un one finds a bifurca- 
tion into three distinct quasi-steady states as introduced 
above. This case will be discussed in detail in Sec. [V] 

The reasons for the occurrence of an SR-like maximum 
of the contrast in terms of the underlying many-particle 
dynamics are illustrated in Fig. [5] To obtain these re- 
sults we have simulated the dynamics generated by the 
Master equation ([7]) using the Monte Carlo wave func- 
tion (MCWF) method (Hll, El H averaging over 100 
quantum trajectories. For a given particle number n, the 
probabilities P to obtain the population imbalance s z 
and the relative phase cj> in a projective measurement are 



So far we have demonstrated a stochastic resonance 
of the contrast for a BEC in a static double-well trap 
with biased particle losses. In the following we will show 
that the system's response to a weak external driving 
also assumes a maximum for a finite dissipation rate - 
an effect which is conceptually closer to the common in- 
terpretation of stochastic resonance. From a mathemat- 
ical viewpoint, however, one can rather relate the un- 
driven case discussed above to the stochastic resonance 
effect in nuclear magnetic resonance (39[. In fact, the 
Bloch equations for the magnetization have constant co- 
efficients in the rotating wave approximation, and should 
thus be compared to the undriven equations of motion 
©. 

Let us consider the response of the system to a weak 
sinusoidal driving of the tunneling rate 

J{t) = J + Ji cos(ujt) (26) 

at the resonance frequency oj = y/ Jq + e 2 , while the am- 
plitude of the driving is small and fixed as Ji/Jn = 10%. 
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(a) 



FIG. 6: Dynamics of the relative population imbalance 
s z (t)/n(t) in a weakly driven double well trap for three differ- 
ent values of the tunneling rate: Jo = 0.5 s _1 (a), Jo = 1.5 s _1 
(b) and Jo = 5s _1 (c). The amplitude of the forced oscilla- 
tions is maximal for intermediate values of Jo as shown in 
part (b). The remaining parameters are T™ 1 — 2s™ 1 , U = 0, 
e = 0, 7 P = 5s™ 1 and Ji/Jo = 10%. Please note the different 
scalings. 



A variation of J can be realized in a quite simple way in 
an optical setup [2l|, [2^, HH , where the tunneling barrier 
between the two wells is given by an optical lattice formed 
by two counter-propagating laser beams. A variation of 
the intensity of the laser beams then directly results in a 
variation of the tunneling rate J. Fig.[5]shows the result- 
ing dynamics for T\ = 0.5 s and three different values of 
Jo and U — 0. After a short transient period, the relative 
population imbalance s z (t)/n(t) oscillates approximately 
sinusoidally. One clearly observes that the response, i.e. 
the amplitude of the forced oscillations, assumes a maxi- 
mum for intermediate values of Jo matching the external 
time scale of the dissipation given by T^ 1 . 

For a detailed quantitative analysis of this stochastic 
resonance effect, we evaluate the amplitude of the oscil- 
lation based on a linear response argument for U = 0. In 
the following, we will use a complex notation for nota- 
tional convenience, while only the real part is physically 
significant. The equations of motion ([8]) are then rewrit- 
ten in matrix form as 



Mie 1 



The matrices Mq and Mi are defined by 



Mo 



/T-T 1 26 

-2e T^ 1 

2J 

V 





-2J 






faTr 



and 



Mr = 




t; 



/aTf 1 Tf 1 / 



0\ 

-2Ji 



o 0/ 



(27) 



(28) 



(29) 



0.15- 




FIG. 7: (Color online) (a) Response (amplitude of the oscil- 
lations of s z (t)/n(t)) of a weakly driven double well trap vs. 
T -1 and Jo calculated within linear response theory, (b) For a 
fixed value of the tunneling rate Jo = 2.5 s" 1 . (c) For a fixed 
value of the dissipation rate T™ 1 = 2s™ 1 . The remaining 
parameters are U = 0, e = 0, 7 P = 5s™ 1 and Ji/Jo = 10%. 



As before we consider the dynamics after all transient 
oscillations have died out, assuming that s(t) as well as 
n(t) decay exponentially at the same rate. However, we 
now also have an oscillating contribution so that we make 
the ansatz 



s(t) = (so+Sie^e"** 
n(t) = (n + me"- t )e™ Kt . 



(30) 



The amplitude of the oscillations, i.e. the system re- 
sponse, is thus directly given by Si/n . Substituting this 
ansatz in the equations of motion and dividing by e _Kt 
yields 



so 



no 

[M + Mie 



(31) 



s \ giwt 

n 



Neglecting the higher order terms ~ e 2»^t m a ii near 
response approximation and dividing Eqn. (|3 1[) in the 
time-dependent and the time-independent parts yields 
the equations 



[-M + (iw - k)1] 



= Mi 



(32) 



and ()15[) . which determine si and ni. The resulting val- 
ues of the system response are shown in Fig. [7J One 
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FIG. 8: (Color online) Dynamics of the coherence s x (t)/n(t) 
(a) and the relative population imbalance s z (t) /n(t) (b) for a 
double well trap with a driven energy bias e for Jo = 2 s _1 and 
Tj -1 = 4s" 1 . (c) Response (amplitude of the oscillations of 
Sx{t)/n(t)) vs. Tj -1 and Jo calculated within linear response 
theory. The remaining parameters are U — 0, ei = Is - , 
7 P = 5s~\ 

observes the characteristic signatures of a stochastic res- 
onance: If one of the two parameters Jo and T\ is fixed, 
the response assumes a maximum for a finite value of the 
remaining parameter as shown in part (b) and (c) of the 
figure. Part (a) shows that this maximum is assumed 
if the external (T^ -1 ) and the internal (Jo) timescale are 
matched similar to the undriven case illustrated in Fig. [2] 
Let us stress that this scenario is again not fundamentally 
altered in the case of weak interactions as numerically 
tested but not shown here. 

A different situation arises if the energy bias is driven 
instead of the tunneling rate J such that 

e(i) = ex cos(ui). (33) 

As above we can evaluate the amplitude of the forced 
oscillations within the linear response theory, however 
with 



Mj = 






-2ei 








2e x 























o 








o 



(34) 



Solving the equations (|3"2")) and (fTTJ)) then yields s Xy = 
Siz = 0. Remarkably, a driving of the energy bias does 




FIG. 9: (a) Time evolution of the purity p and the contrast 
a for J = U = 10 s" 1 , e = 0, Ti = 0.5 s. (b) Time evolu- 
tion without interactions (U — 0) and (c) without dissipation 
(1/Ti = I/T2 = 0) for comparison. The occasional revivals 
are artifacts of the small particle number. The initial state 
is a pure BEC with s z — n/2 and n(0) = 100 particles. The 
results of a MCWF simulation averaged over 100 runs are 
plotted as a thin solid line in (a) and (c) , while the mean-field 
results are plotted as a thick line in (a) and (b). Note that 
the mean- field approximation is exact in case (b), whereas it 
breaks down in case (c) and is thus not shown (cf. [291. 1321]). 

not affect the population imbalance in leading order. 
Only the first component of the Bloch vector s x , and 
thus also the contrast a is strongly affected. 

This is illustrated in Fig. [5] (a) and (b) where the rel- 
ative population imbalance s z (t)/n(t) and the first com- 
ponent of the Bloch vector s x (t)/n(t) are plotted for 
J = 2s" 1 , Tf 1 = 4s" 1 and e x = Is -1 . The coher- 
ence oscillates strongly at the fundamental frequency w, 
while the population imbalance oscillates only with a tiny 
amplitude at the second harmonic frequency 2u>. The os- 
cillation amplitude of the coherence then again shows the 
familiar SR-like dependence on the parameters Jo and T\ 
as illustrated in Fig. [8] (c). 

V. DISSIPATION INDUCED COHERENCE IN A 
STRONGLY-INTERACTING BEC 

Let us finally discuss the case of strong interactions, 
which is experimentally most relevant and theoretically 
most profound. An example for the dynamics of a 
strongly-interacting BEC is shown in Fig. [5] (a) for an 
initially pure BEC with s z = n/2, calculated both with 
the MCWF method and within the mean-field approxi- 
mation (jSJ) . One observes that the purity p and the con- 
trast a first drop rapidly due to the phase noise and, 
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FIG. 10: (Color online) Mean-field dynamics without interac- 
tions and dissipation (a), with interactions Un = 40s -1 (b) 
and with interactions and dissipation j a = 10s _1 (c). The re- 
maining parameters are J = 10s _1 and e = 0. To increase the 
visibility we have plotted the rescaled Bloch vector s/n and we 
have artificially fixed the particle number so that n = const. 



more importantly, due to the interactions. This is an 
effect well-known from the non-dissipative system and 
can be attributed to a dynamical instability which also 
leads to the breakdown of the mean-field approximation 
[29l . [30l . l32l . Hi]]. However, a surprising effect is found 
at intermediate times: The purity p is restored almost 
completely and the contrast a is slightly increasing. 

Most interestingly, the observed values of the purity 
and the coherence are much larger than in the cases where 
one of the two effects interactions and dissipation - is 
missing. The time evolution for these two cases are also 
shown in Fig. [5] In the case of no interactions both purity 
and coherence rapidly drop to values of almost zero and 
do not revive. This case has been discussed in detail 
in Sec. IIIII In the interacting case without dissipation 
one observes regular revivals, which are artifacts of the 
small particle number in the simulation and become less 
pronounced with increasing particle number. Apart from 
these occasional revivals, however, the purity and the 
coherence relax to values which are much smaller than in 
the interacting and dissipative case. 

The surprising re-purification of a strongly-interacting 
BEC by particle dissipation can be understood within a 
semiclassical phase space picture. In order to visualize 
the effects of particle loss, we have plotted the 'classical' 
phase space structure generated by the Bloch equations 
((HI) for 7 P = in Fig. [TO] without interactions and dissi- 
pation (a), with interactions (b) and with both (c). For 
illustrative purposes, we have plotted the rescaled Bloch 
vector s/n and have artificially fixed the particle num- 
ber so that n — const. Since we are interested only in 
the short-time dynamics of the Bloch vector and not in 




0123401 234 
t [s] t [s] 

FIG. 11: Time evolution of the purity p and the contrast a for 
J = U = 10s" 1 , e = and 1/Ti = 0.5s" 1 (a), 1/Ti = 1.5s" 1 
(b) and 1/Ti = 2.5s" 1 (c). The initial state is a pure BEC 
with s z = n/2 and n(0) = 100 particles. The results of a 
MCWF simulation averaged over 100 runs are plotted as a 
thin solid line, while the mean-field results are plotted as a 
thick line. 



the decay of the particle number on longer time scales, 
this is an appropriate treatment. Moreover, in the quan- 
tum jump picture this approximation corresponds to the 
periods of constant particle number between two loss pro- 
cesses [HiHiiHii. 

Parts (a) and (b) of the figure show the phase space 
structure without dissipation and Un = and Un = 4 J, 
respectively. One observes the familiar self-trapping bi- 
furcation of the fixed points for Un > 2 J [28|, |29| . The 
phase space structure is significantly altered in the pres- 
ence of particle loss as shown in part (c). The most im- 
portant consequence is the occurrence of an attractive 
and a repulsive fixed point instead of the elliptic fixed 
points in the dissipation- free case [38| . 

In the course of the time the system will thus relax 
to the attractive stationary state illustrated Fig. [TO] (c). 
A many-particle quantum state can now be represented 
by a quasi-distribution function on this classical phase 
space, for instance the Husimi Q-function [3l], [32[. In 
this picture, a pure BEC is represented by a maximally 
localized distribution function and the loss of purity cor- 
responds to a broadening or distortion of the Q-function. 
The existence of an attractive fixed point clearly leads 
to the contraction of a phase space distribution function 
and thus to a re-purification of the many-particle quan- 
tum state as observed in Fig. [5] (a). 

However, this nonlinear stationary state exists only as 
long as the particle number exceeds a critical value given 
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(a) 1 




FIG. 12: Purity p (a) and contrast a (b) after t = 2 s as a 
function of the dissipation rate 1/Ti for different values of 
the interaction strength Un calculated within the mean-field 
approximation. The remaining parameters are chosen as in 
Fig. M (a). 

by (cf. [U) 

U 2 n 2 > 4 J 2 - / 2 Tf 2 . (35) 

As particles are slowly lost from the trap, the particle 
number eventually falls below the critical value. For this 
reason the attractive fixed point vanishes and the purity 
drops to the values expected for the linear case U = 0. 
Since the attractive fixed point tends towards the equator 
maximizing s^/js), the contrast assumes a maximum just 
before the disappearance of the attractive fixed point, 
while the purity is still large. In Fig. [5] (a) this happens 
after approx. 2.5 s. 

The surprising effect of the re-purification of a BEC 
is extremely robust - it is present as long as the condi- 
tion (|35[) is satisfied. A variation of the system param- 
eters does not destroy or significantly weaken the effect, 
it only changes the time scales of this relaxation process. 
Fig. [IT] compares the time evolution of the purity and 
the contrast for three different values of the particle loss 
rate Tf \ With increasing losses, the nonlinear station- 
ary state is reached much faster, but is also lost earlier. 
One can thus maximize the purity or the contrast at a 
given point of time by engineering the loss rate. This ef- 
fect is further illustrated in Fig. 1121 where the purity and 
the contrast after 2 seconds of propagation are shown in 



dependence of the loss rate T l 1 . Both curves assume a 
maximum for a certain finite value of T± . 

VI. CONCLUSION AND OUTLOOK 

In summary, we have shown that the coherence prop- 
erties of a weakly and, in particular, also of a strongly in- 
teracting Bose-Einstein condensate in a double-well trap 
can be controlled by engineering the system's parameters 
and dissipation simultaneously. Surprisingly, dissipation 
can be used to stimulate coherence in the system, rather 
than - as may be expected - solely reduce it. 

In the weakly interacting case, the contrast of the 
quasi-steady state of the system assumes a maximum for 
a finite value of the tunneling and the dissipation rate. 
This stochastic resonance effect is robust against param- 
eter variations. A Monte Carlo wave function simulation 
of the full many-body quantum dynamics shows a good 
agreement to the mean-field description and provides a 
microscopic explanation of the observed effect. More- 
over, a similar effect can be observed in the case where 
either the tunneling or the energy bias is driven, which 
is conceptually even closer to the common interpretation 
of stochastic resonance. 

In the last section, we have studied the effects of dissi- 
pation on the strongly interacting system. An important 
conclusion is that the interplay of interactions and dissi- 
pation can drive the system to a state of maximum co- 
herence, while both processes alone usually lead to a loss 
of coherence. We show that this effect can be understood 
from the appearance of an attractive fixed point in the 
mean-field dynamics reflecting the metastable behaviour 
of the many-particle system. 

Since the double-well BEC is nowadays routinely re- 
alized with nearly perfect control on atom-atom interac- 
tions and external potentials [2l|, [2^, [Hj], we hope for 
an experimental verification of the predicted stochastic 
resonance effect. An interesting perspective is to lift our 
results to extended dissipative setups, as e.g. studied in 
[551 . [56[. Besides the general idea of controlling many- 
body dynamics [57], one may also investigate the possi- 
bility of dynamically engineering entanglement in similar 
systems, as to some extend possible in state-of-the-art 
experiments (58|. 
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